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Abstract 

The form factors and the coupling constant of the BsB*K and BgBK* vertices are calculated 
using the QCD sum rules method. Three point correlation functions are computed considering 
both the heavy and light mesons off-shell in each vertex, from which, after an extrapolation of the 
QCDSR results at the pole of the off-shell mesons, we obtain the coupling constant of the vertex. 
The form factors obtained have different behaviors but their simultaneous extrapolation reach the 
same value of the coupling constant gBsB*K = 8.41 ± 1.23 and gSsBK* = 3.3 ± 0.5. We compare 
our result with other theoretical estimates and compute the uncertainties of the method. 
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I. INTRODUCTION 


In the recent years, many new charmonium and bottomonium states have been observed 
at the B-factories. As an example, in the bottomonium sector, the Belle Collaboration 
reported the observation of two charged narrow structures in the 7r^T(nS') {n = 1, 2, 3) and 
'K^hhijnP) {m = 1,2) mass spectra of the T(5S') —)■ T(nS')7r=‘= and T(5S') —)■ hi,{mP)n^ 
decay processes [T]. These narrow structures were called Zfe(10610) and Zb(10650). As 
pointed out by the Belle Collaboration, the proximity of the BB* and B*B* thresholds 
and the ^^(10610) and Zf,(10650) masses suggests that these states could be interpreted as 
weakly bound BB* and B*B* states. In particular, using the one-boson exchange model 
and considering S'-wave and H-wave mixing, the authors of Ref. |2] were able to explain 
both, Zf,(10610) and Zft(10650), as BB* and B*B* molecular states. The main ingredients 
in the one-boson exchange model are the effective Lagrangians, that describe the strong 
interactions between the heavy and light mesons. These Lagrangians are characterized by 
the strong coupling constants in the considered vertices which, in general, are not known. 
These heavy-heavy-light mesons coupling constants are fundamental objects, since they can 
provide essential information on the low energy behavior of the QCD. Depending on their 
numerical values, a particular molecular state may or may not be bound. Therefore, it is 
really important to have reliable ways to extract these values based on QCD calculations. 
However, such low-energy hadron interaction lie in a region which is very far away from 
the perturbative regime. Therefore, we need some non-perturbative approaches, such as the 
QCD sum rules (QCDSR) |3H5], to calculate the form factors and coupling constants of these 
vertices. There are already some QCDSR calculations for the heavy-heavy-light vertices like 
the B*B7t [6], BsoBK [7|, B*BK* [8], 0, BsBK*, B*BKi [TO] and B*BK [H]. 

In the charm sector, various vertices were evaluated with this approach and the results are 
systematized in [TO]. Here we calculate the form factor and the coupling constant at the 
BgB*K and BgBK* vertices in the framework of three-point QCDSR. More specihcally, we 
evaluate the gBsB*K{Q‘^) and gBsBK*{Q‘^) form factors in three different ways, considering, 
one by one, each one of the mesons in the vertex to be off-shell. From these form factors, 
we extract the gBsB*K and the gB^BK* coupling constants. 
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II. THE QCD SUM RULE FOR THE BsB*K AND B^BK* VERTICES 


To perform the QCDSR calculation and obtain the form factors and coupling constants 
of the vertices BsB*K and BgBK*, we follow our previous works, as Ref. [12]. The starting 
point is the three-point correlation function given by : 


n 


!S) (p, p) = I (0| T ' (0)} |0) 


,kM 


t/. 


o-n-y 


( 1 ) 


for the Bs meson off-shell, and: 

= J d^xdS | 0 ) ( 2 ) 

for the K or K* meson off-shell. In Eqs. Q and Q, g = p — p' is the momentum of the off- 
shell meson and p and p' are the momentum of other ones. The currents are the currents 
associated with each meson in the vertex and contain the quantum information about the 
state. The correlation functions in Eqs. Q and ([^ allow to obtain two different form factors 
corresponding to the same vertex. In this way, the vertex is tested by two different mesons, 
the heavier and the lighter ones in the corresponding vertex. The calculation of these two 
correlation functions allows to reduce the uncertainties of the evaluation of the coupling 
constant of the vertex na. 

Equations Q and ([^ contain different numbers of Lorentz structures, and for each struc¬ 
ture, we can write a different sum rule. In principle all the structures would give the same 
result. However, due to different approximations each structure can lead to different results. 
Therefore, one has to choose the structures less sensitive to the different approximations. 
To obtain the sum rule, these functions are calculated in two different ways: using quarks 
degrees of freedom - the QCD side; and using hadronic degrees of freedom ~ the phenomeno¬ 
logical side. In the QCD side, the correlators are evaluated using Wilson’s operator product 
expansion (OPE). The duality principle allows us to obtain an interval in which both repre¬ 
sentations are equivalent. Therefore, in this region, we can obtain the QCD sum rule from 
where the form factors are evaluated. To improve the matching between the two sides, we 
perform a Borel transformation to both QCD and phenomenological sides. 
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A. The QCD side 


The QCD side is obtained using the following meson currents for the BsB*K vertex: 

jT (^) = 

jgS-(O) = i675s, 

Juiy) = (3) 

and the following ones for the BgBK* vertex: 

ji{x) = iqibh, 

jf-(O) = z575s, (4) 

Here, q, s and b are the light, strange and bottom quark helds respectively. Each one of 
these currents has the quantum numbers of the associated meson. In the case of K off-shell 
meson, we use, as usual in QCDSR, the pseudo scalar current for it, see refs. D*Dti [B], 
D^DK [13] and B^BK HU. The general expression for the vertices has different structures, 
which can be written in terms of a double dispersion relation over the virtualities and 
holding = —q^ hxed: 




47^2 


ds 


du 


p(s,m,Q2) 


(s — p^){u — p'^) 


(5) 


where the spectral density p{s,u,Q‘^) can be obtained from the Cutkosky’s rules. 

The invariant amplitudes receive contributions from all terms in the OPE. In the case of 
form factors, the main contribution in the OPE is the perturbative term, which is represented 
in Fig. for the two cases that we are considering, the Bg^Bg) and K(K*) meson off-shell 
for the BgB*K (BgBK*) vertex. 

In order to obtain the form factor, we have to choose one of the different structures ap¬ 
pearing in Eqs.([^ and ([^. As commented above, different structures can lead to different 
results. Therefore, one has to choose the structure less sensitive to higher dimension con¬ 
densates, that provide a better stability as a function of the Borel mass, and that have a 
larger pole contribution, when compared with the continuum contribution. This is consid¬ 
ered a “good” structure. If there is more than one “good” structure, the others can also be 
considered to estimate the uncertainties of the method. 
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FIG. 1: Perturbative diagrams for the Bs{Bs) off-shell meson (left) and for K{K*) off-shell meson 
(right), for BsB*K vertex {BgBK* vertex). 

For BsB*K form factor and in the case Bg off-shell meson, we choose the structure, 
because it satishes the criteria above. For K off-shell meson, we can work with both p^ and 
p'j^ structures In this case, we are going to work with the p), structure while the other one 
will be used for the estimate of the uncertainties. The corresponding perturbative spectral 
densities, which enter in Eq. ([^, are : 


u, Q^) = -^ [{2mb - 2ms)E - 2mbB] , 


for the p'^p'j, structure of the Bg off-shell case, and 




271^/X L 


27iy/X 

e, and 

A{p ■ p — 2k ■ p — mbmg + ml) + 27i{ml — k ■ p') 


( 6 ) 


(7) 


for the py structure in the K off-shell case. In both cases, the quark condensate contributions 
are neglected after the Borel transform. 

For BgBK* vertex, we use the p^ structure for Bg off-shell meson and for K* off-shell 
meson, we have p^ and p'^ structures, both giving excellent sum rules. Again we show the 
results for one structure and the other is used to estimate the uncertainties. The perturbative 
contribution to the spectral density, when the Bg meson is off-shell is : 




27ryA 


[A{p ■ p' — mbmg — 2p ■ k) — p ■ k] 


( 8 ) 


for the p^ structure. In this case, the quark condensate, (qq), contribution to the same 
structure is: 


n = 


mg{qq) 


(p2 — ml){p'‘^ — ml) 

For K* off-shell case, the spectral density, for both structures, is given by: 
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( 9 ) 


p^^*\s,u,Q‘^) = 


2tt\/~X 


{p^ A{p ■ p - mbmg - ml) + ml - k ■ p - mbmg) 
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+ p'fM B{p ■ p'+ rahTiis — ml) — k ■ p + ml }. (10) 

In Eqs. (|^ to (10), we have defined A = A(s, u, t) = — 2st — 2su — 2tu, s = 

u = p'^, t = —Q"^ and A, B and E are functions of {s,u,t), given by: 


A = 271 


ko ^ 


s \pA \/s 


Po 


\k\ 


B = 271==cos9 ; 


( 11 ) 


E = (3cos0- l) , 

|p-| ^ ^ 

where 


( 12 ) 
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(13) 


Finally, the OPE side, is calculated using Eq. (|^ with the limits in the integration given 
by: Smin = ijnb + ms)'^ and Umin = t — ml for Bg off-shell and Smin = f^l — ml and 
Umin = t + ml — ml for K off-shell, for BsB*K vertex. And Smin = {'rrib)'^ and Umin = t + ml 
for Bg and K* off-shell for the BgBK* vertex. 


B. The phenomenological side 

The three-point functions from Eqs. 0 and (|^ , when written in terms of hadron masses, 
decay constants and form factors, give the phenomenological side of the sum rule. 


1. For the BgB*K vertex: 

The meson decay constants fx, fss and fs* are defined by the following matrix elements: 


(ob7|A:(p)) = ifxPu, 

(14) 

(0|jf|i?*(p'))=rnB*/B*6;(p'), 

(15) 

mb + mg 

(16) 
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and the vertex function is defined by: 


{K{q)\B-(-p')B.(p)) = -3^B.l.K^^'{p')(2p-p'),, 
which is extracted from the effective Lagrangian HU: 


(17) 


C 


BsB*K — WBsB*K 


B'^(B,d^K - d^B.K) - B’HBsSbK - d„B.K) 


(18) 


Saturating the correlation function with Bg, B* and K intermediate states we arrive at 

-fKfBsfBmB*m\Q^ZB*K{(f) 


— _ 

{irts + mb){p'‘^ — m\){q2 — 'rn%J{p‘^ — m^*) 


X 




1 _ vx 


m 


B* 


flB V 


+ “continuum’ 


for an off-shell Bg. Using the matrix element of K meson equal to 

TRg 

we arrive at an expression for an off-shell K: 

^ -fB*fKfBBmB*m%m\g^BsB*K{<f) 

^ {mb -|- mg)mg{p‘^ — m\){pf'^ + m\,){q^ — mj^) 


-‘^Pu + P u 1 + 






-|- “continuum’ 


(19) 


( 20 ) 


( 21 ) 


2. For the BgBK* vertex: 


In this case, the effective Lagrangian is |14j : 


C 


BbBK* — 'I'QBsBK* 


K*^{Bd^Bg - Bgd^B) + K-^^{Bg^^B - Bd^Bg 

from where we can extract the vertex element, which is given by: 

(K-(q)\B,{-p')B{p))=ig'-^;2l<.r‘‘{q)(p + p‘)B, 


( 22 ) 


(23) 


and the matrix elements which introduce the meson decay constants fx*, and fs are: 


= fK*e*^{q)mK*, 


(24) 


{B{p)\ji\0) 


m 




mt 


(25) 
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and 
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Bs 


m^\Bs{p^)) = fs. ^ , 

TTlb + rUg 

After some algebra we arrive at the following expression: 


(26) 


ll(Bs) ^ 
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{ml + msmb){p^ — m%){pf‘^ — m\t){q^ — m 
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+ “continuum”, 


(27) 


when Bg is off-shell. 


For K* off-shell we arrive at: 




fBfK-fBsmK*rnl m 
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^9fsBK*{(f) 
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+ “continuum’ 


(28) 


III. THE SUM RULE 

The snm rnle is obtained after performing a donble Borel transform {BE), = —p^ —)■ 

and P'^ = ——)■ to both the phenomenological and OPE sides: 


BE 


-nOPE{I) 


(M, M') = BE 


Y'phen^I) 


(M,M'), 


(29) 


where M and M' are the Borel masses and I is the off-shell meson. 

In order to eliminate the continnnm contribntion in the phenomenological side, instead 
of doing the integrals in Eq. ([^ np to cx), we do the integrals np to the continnnm threshold 
parameters sq and uq. The threshold parameters are dehned as sq = and uq = 

(mo -|- Ao)^, where Aj and Ao are usually taken as 0.5 MeV, and m* and mo are the masses 
of the incoming and outgoing mesons respectively. 

In Eqs. (19), (21), (27), (28), g^Bls^RiQ'^) gBsBK*iQ‘^) the form factors when the 
I meson is off-shell. As in our previous works, we define the coupling constant as the value 
of the form factor, at = —ml where mj is the mass of the off-shell meson. 




















IV. RESULTS AND DISCUSSION 


Table shows the value of the hadronic parameters used in the present calculation. We 
have used the experimental value for fx of Ref. [13] , for fs* and from Ref. [TB] and for 
fx* of Ref. [H] . 


TABLE I: Hadronic parameters used. 


K 

K* 

Bg 

B* 

B 

m (GeV) 0.49 

0.89 + 5 

5.40 

5.20 

5.28 

/ (MeV) 

220 + 5 

208|io 

250^29 

191 + 0.87 


We neglect the light quark mass {niq = 0.0 MeV). The strange and bottom quark masses 
were taken from the Particle Data Group (PDG) and have the values mg = 104+26 —34 MeV 
and mfc = 4.20 + 0.17 — 0.07 GeV respectively. In the next subsections, we show the two 
form factors used to extract the coupling constants of each vertex. 

A. BsB*K vertex 




(b) 


FIG. 2: (a) The pole-continuum contributions for the BsB*K sum rule for the Bg off-shell and (b) 
9i^sB*KiQ‘^ = 1 GeV^) stability for different values of the continuum thresholds. 


In the case of a Bg off-shell meson, we work with the structure. In Fig. j^a) we 
show the contribution of the pole versus the continuum contribution for the sum rule and in 
Fig. I^b) the stability of the form factor as a function o Borel mass, for = 1 GeV^ and 


9 













three different values for the continuum thresholds. We use the usual relation between the 
Borel masses M' and M |T21: M'^ = — 

From Fig. |^b) we clearly see a window of stability for > 10 GeV^, and from Fig. [^a), 
we can see that the pole contribution is bigger than the continuum contribution for < 
19 GeV^. Therefore, there is a Borel window where the sum rule can be used to extract the 
form factor. 


Pole 

Continuum 


M^(GeV^) 

(a) 



(b) 


FIG. 3: a) The pole and continuum contributions for BsB*K sum rule with the K off-shell and b) 
the ~ ^ GeV^, M^) stability for different values of the continuum thresholds 


In the case of a iF meson off-shell, we have two structures in Eq. (21) that can be 
used: p'^ and Both structures give good sum rule results, that means a good pole- 
continuum contribution and good stability. We show in Fig. [^a) and (b), repectively, the 
pole-continuum contribution and the stability only for the p^, structure. For the p'^, we 
obtain a very similar result and we use it to evaluate the uncertainties. 

From Fig. we find a Borel window 10 GeV^ < < 14 GeV^ where the sum rule can 

be used to extract he form factor. 

In Fig. [^we show the QGDSR results for these two form factors, represented by squares 
and triangles for the cases of the Bs and K off-shell respectively. 


B. BgBK* vertex 

For the case of the Bg meson off-shell, we use the structure p^. In Fig. |^a) and (b) we 
show the pole-continuum contribution contribution and the Borel mass stability respectively. 
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FIG. 4: a) The pole-continuum contributions for BgBK* and in b) ~ lGeV^,M^) 

stability for different thresholds. 

Once more, we see that we get a Borel window lOGeV^ < < 12.5 GeV^ where these 

both conditions are satisfyed and where we can use the sum rule to extract the form factor. 
Also, for the Bg meson off-shell case, we have other structure, p'^, to work with. In this case, 
the result includes the light quarks condensate which do not modihed substantially the sum 
rule. The p'^ structure will be used to evaluate the uncertainties. 

In the case of the K* off-shell meson, we show the hgures of stability and pole-continuum 
contribution for p'^ structure, in Figj^ In the structure we get similar results, that means 
a good pole-continuum contribution and good stability. 




FIG. 5: a) pole and continuum contributions for BgBK* sum rule for K* off-shell and b) the 
g^BsBK*iQ‘^ = 2 GeV^, M^) stability for different values of the threshold parameters 

The sum rule result for both form factors are shown in Figthrough triangles and squares 
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for Bs and K* off-shell respectively. 


V. COUPLING CONSTANT 


The coupling constant in the vertex, is obtained when the form factor is extrapolated 
to = —mf, where m/ is the mass of the off-shell meson. In order to minimize the 
uncertainties, we work with two form factors in each vertex, one with the heavier meson in 
the vertex off-shell and the other with the lighter meson off-shell. Both form factors should 
give the same value for the coupling constant. 

In Table [11} we show, for each vertex, the function that hts the QCD sum rule results 
for the form factors, which is then extrapolated to = —mj, to determine the coupling 
constant value. From Table we see that the two form factors in the same vertex, give 
similar results for the coupling constant. 


TABLE II: Form factors and coupling constants obtained by extrapolation 
Form Factor Coupling Constant 


= 3.18exp <3^/32-98 qb^b-k = 7.55 

(Q') = 8-25 exp-(«^/6-68)^ gB^B*K = 8.00 

9BsBK*(Q'^) = 58*86+Q^ 9BsBK* = 2.72 

9^BsBK*iQ‘^) = 2.49exp-'3V2.63 qb.bk* = 3.37 


In Figs. [^and[^ we show the parametrizations given in Table [IT| for the QCDSR results 
for the form factors gBsB*K{Q^) and gBsBK*{Q‘^) respectively. The squares and triangles give 
the QCDSR results for the heavier and lighter off-shell mesons in the vertices respectively. 
We also show in these hgures the values obtained for the coupling constant. 

VI. UNCERTAINTIES OF THE QCDSR METHOD 

Looking at Fig. |^and[^ we can see a error bar at the endpoint of the curves. To evaluate 
these error bars, we compute the sum rule taking into account the errors in the masses, decay 
constants, continuum threshold parameters and also we study a variation on the Borel mass 
in the window of Q^. 
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FIG. 6: The gBsB*K{Q‘^) form factor, extrapolation function for K and Bs meson off-shell. 



FIG. 7: The gBsBK*{Q‘^) form factor, extrapolation function for K* and Bg meson off-shell. 


Our procedure is that in each computation, all the parameters are kept hxed and only 


one changes within its intrinsic error. In Table |III[ we show the intrinsic uncertainty for 
each of the sum rule parameters. 

We also consider the other good sum rules, for each off-shell meson in each vertex and 
we also put the third meson off-shell to better estimate the errors. 
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TABLE III: Percentage deviation of the coupling constant related to each parameter for both 
vertices with the meson I off-shell. 



Deviation % 

SBsB*K UBs 

BK* 

Parameters 

I = Bg 

I = K 1 

= Bg 

I = K* 

Jk* = 220 ± 5(MeV) 

~ 

- 

1.95 

1.94 

fB = 191 ± 8.77(MeV) 

- 

- 

3.91 

3.88 

Jk = 159.8 ± 1.4 ± 0.44(MeV) 

1.04 

1.04 

- 

- 

fB* = 208 ± 10 ± 29 (MeV) 

16.29 

16.30 

- 

- 

/s, = 250 ± 10 ± 35 (MeV) 

12.06 

12.07 

3.61 

3.51 

mb = 4.20 -7 0.17 - 0.07 (GeV) 

15.81 

23.39 

9.83 

28.42 

mg = 104 26 - 34 (MeV) 

1.23 

24.71 

7.21 

19.78 

m2 ± 10% (GeV) 

4.29 

1.85 

6.45 

18.57 

As ± 0.1 e An ± O.l(GeV) 

13.51 

4.90 

9.32 

12.35 


A. Other mesons of shell and other structures 

For the BsB*K vertex, we calculated the B* off-shell form factor. Since the B* meson 
has a small mass difference when compared with the Bg meson mass, we expect to obtain 
a similar result to the one obtained for Bg off-shell. Computing the sum rule, we observe 
that there are two structures: the and p^p^. The p^p'^ structure does not give a good 
stability in the Borel mass, therefore it is not considered. The other structure, p'^p'i^, gives a 
good stability and the pole is larger than the continuum contribution. 

To extrapolate the QCDSR results, we use the £t 

= (30) 

for the form factor, and the resulting value for the coupling constant is: gBsB*K = 8.64. This 
result is very similar to the one showed in Table as expected. 

Also for this vertex, in the case of the K meson off-shell, we have other structure, p'^, 
which gives a good sum rule, very similar to the p^ used before. In this case the form factor 
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is extrapolated by 



( 31 ) 


and the resulting value for the coupling constant is: gBsB*K = 9.45, in a very good agreement 
with the results in Table m 

For BgBK* , we perform the sum rule for the B off-shell meson. Again, we expect to 
obtain a similar result to the one obtained for Bg off-shell, since this meson has a small mass 
difference when compared with the B meson mass. 

We get sum rules for the two structures, and p'^. The structure does not give good 
stability in the Borel mass. Therefore, it is not considered. The other structure, p'^, gives 
good stability with the pole being bigger than the continuum contribution. The QCDSR 
results for the form factor can be extrapolated by 



(32) 


and the coupling constant is 


QBsBK* — 2.77. 


This is a very similar result to the one obtained in Table 

VII. CONCLUSION 

We have used the three-point QCD sum rules to study the form factors in the vertices 
BgB*K and BgBK*. In each case, we have considered two different sum rules, for two dif¬ 
ferent mesons off-shell. We have studied the Borel stability and pole dominance in each case 
and have determined the Borel window where the sum rules can be used. To get the coupling 
constant in each vertex, we have used the extrapolation method developed in previous works 
nzi. This extrapolation method has a systematic error which comes from the choice of the 
analytic form of the extrapolating functions. We consider only monopole, exponential and 
Gaussian parametrization: 5 'mm 2M3(Q^) = Aexp“‘3AB^ 9mm2M3{Q‘^) = Aexp'^'^VB)^ 

There are no physical reason for using these forms. However, they are the usual forms 
used by experimentalists and also they have only two parameters, A and B, that present 
some regularity, in all our form factors. For instance, when the heavy meson is off-shell, 
the form factor is harder as a monopolar function, and the B cut-off parameter is a big 
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number. By the other hand, when the light meson is off-shell, the parametrization is usually 
exponential or Gaussian and have a smaller B cut-off parameter. These cut-offs are showed 
in Table IB 

We have also performed a very extensive study of the uncertainties to estimate the errors 
in the coupling constants. After considering all good sum rules and including the uncertain¬ 
ties study, we obtain the coupling constant equal to: 

9b.b*k = 8.41 ± 1.23 (33) 

and 

Qb^bk* = 3.3 ± 0.5 (34) 


cut-ofF parameters 



BsB*K 

BsBK* 

Off-shell meson 

A B 

A B 

Bs 

3.18 32.98 

80.88 58.86 

B* 

1.21 13.65 

- 

K 

8.25 6.68 

- 

B 

- 

104.71 66.97 

K* 

- 

2.49 2.63 


TABLE IV: Cut-off parameters for both vertices. 

We can compare our results with other theoretical predictions using arguments of heavy 
hadron chiral perturbation theory (HHChPT), where the couplings for the bottom-light 
vertex Qb^bk are related to the charm-light vertex gosO^K through the relation |THl HH]: 

gBsB*K = gD%DK — , (35) 

rriD 

where = 5.2 and mn = 1.8693 are the experimental masses. For Qd^d^k we can use our 
previous QCDSR result [13] of gD*DK = 3.02. Therefore, we obtain gBsB*K = 8.39, which 
is in excelent agreement with our result. This result shows that the relation with charm- 
light vertex can give a good estimate of the couplings in the bottom-light vertex. The 
uncertainties are about 20%, that is in complete agreement to the technique of QCDSR. 
The greatest source of uncertainties found in this work is due to the mass of the bottom 
quark. 
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